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We present a very general argument that the analogue of a heavy pentaquark (a state with the 
quantum numbers of a baryon combined with an additional light quark and a heavy antiquark Q) 
must exist as a particle stable under strong interactions in the combined heavy quark and large 
Nc limits of QCD. Moreover, in the combined limit these heavy pentaquark states fill multiplets 
of SU(4) X 0(8) X SU(2). We explore the question of whether corrections in the combined 1/Nc 
and 1/mQ expansions are sufficiently small to maintain this qualitative result. Since no model- 
independent way is known to answer this question, we use a class of "realistic" hadronic models 
in which a pentaquark can be formed via nucleon-heavy meson binding through a pion-exchange 
potential. These models have the virtue that they necessarily yield the correct behavior in the 
combined limit, and the long-distance parts of the interactions are model independent. If the long- 
^ , distance attraction in these models were to predict bound states in a robust way (i.e., largely 

insensitive to the details of the short-range interaction), then one could safely conclude that heavy 
pentaquarks do exist. However, in practice the binding does depend very strongly on the details of 
^SJ , the short-distance physics, suggesting that the real world is not sufficiently near the combined large 

Nc, mq limit to use it as a reliable guide. Whether stable heavy pentaquarks exist remains an open 
fT^ ' question. 
> 

Q\ , PACS numbers: 12.39.Mk, ll.15.Pg, 12.39.Hg, 12.39.Pn 
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' The existence of pentaquarks remains a vexing unresolved experimental question. Ten groups performing a variety 
of experiments P,!^^] have reported the appearance of the pentaquark state now called 8+, a resonance with baryon 
O • number strangeness -t-1, and a mass in the vicinity of 1540 MeV. However, these experiments were all performed 
I with relatively limited statistics and significant cuts, raising the possibility that the reported resonance is due to 
nothing more than statistical fluctuations. One ground for skepticism arises from a series of experiments that did 
(— I not find a 8+ resonance ^|^. Of course, it is unclear whether some of the experiments with negative results should 
' be sensitive to such an observation, since there is no reliable theoretical framework for predicting the 0"*" production 
rate. The 0+ width generates another source of doubt: r(0+) must be exceedingly narrow (in the range of 1-2 MeV 
or smaller), or it would have been detected long ago 6j, and to many it strains credulity that such a narrow state 



exists in this kinematic range. 

One common thread in these early reports of detection (or non-detection) of the 0^ is the dependence of the 
experimental analysis upon revisited old data, and the appearance of the signal only after the imposition of various 
cuts. Given the limited size of these old data sets, all of the studies yielded spectra with very limited statistics, creating 
the possibility of narrow peaks due to statistical fluctuations. The need for high-statistics experiments became very 
clear. Special-purpose experiments designed to look for pentaquarks with high statistics have been performed at 
Jefferson Lab; the CLAS Collaboration has analyzed the high-statistics data from photons on both a proton target Q 
and a deuterium target j8j, and finds no evidence for a 0+ peak. While these experiments alone do not rule out the 
0+, they show that at least two of the previous claims of evidence for the state, the SAPHIR 77? result and the 
CLAS 7c? result 01, were indeed statistical fluctuations. The prospect that other claims of evidence for the O"*" may 
also evaporate weighs heavily on the field. The initial observation j2j of S pentaquark states appears to be headed 
for a similar fate |5j- 

The theoretical landscape for pentaquarks has been just as murky. A paper by Diakonov, Petrov, and Polyakov 
was seminal in focusing attention on the pentaquark, in that it predicted a narrow state at almost exactly the mass 
where the 0^ was later reported. However, this paper is based upon an approximation later shown to be inconsistent 
with the large assumptions implicit in the model After the experimental claims of pentaquarks appeared, a 
vast literature of models for the 0+ followed. In all of these models the existence of the 0+ depends upon ad hoc 
assumptions; thus they cannot be used reliably to predict the existence of the state, and accordingly are not reviewed 
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here. Ultimately one may hope for lattice QCD eventually to resolve the theoretical question of whether the state 
exists. However, current lattice simulations for both heavy and light pentaquarks |l2l |. while always improving, still 
remain inconclusive. 

Given this morass, it is sensible to ask whether one can find a regime in which the question of t he p entaquark's 
existence is more tractable. It has been noted previously in the context of various models 0, Q| that heavy 
pentaquarks, states in which the s quark in 0+ is replaced by a c or a 6 quark, is more likely to be bound than the 
s type. The principal purpose of this paper is to explore the possible existence of heavy pentaquarks. We show in a 
particular limit of QCD, the combined large Nc and heavy quark limits, that heavy pentaquarks must exist, that they 
are stable under strong interactions, and that they fall into multiplets of SU(4) x 0(8) x SU(2). Here, the SU(4) is the 
large Nc spin- flavor symmetry of the light u and d quarks llfil Il7l |. the 0(8) is a dynamical symmetry associated 
with collective vibrations of the heavy antiquark Q (mass mg) relative to the remainder of the system (18> |. and the 
SU(2) is the symmetry of separate rotations of the Q spin. We then explore the critical question of whether 1/Nc 
and l/mg corrections are sufficiently small for this qualitative result to survive in the physical world. There are no 
known analytic methods starting directly from QCD to answer this last question; thus, we investigate the question in 
the context of models. 

We employ models that treat the heavy pentaquark as a bound state of a heavy meson and a nucleon interacting 
via pion exchange. Although similar models have been considered previously |l9| . the present work expands on them 
and is done in the context of the combined heavy quark and large Nc limits. Such models have two principal virtues: 
First, as we show below, the combined large Nc and mg limit mandates the existence of bound pentaquarks. Indeed, 
our demonstration is based on the fact that QCD in the combined limit can be reduced to a model of this form. 
Second, the long-distance behavior of the model is well known empirically (up to experimental uncertainties in the 
pion- heavy meson coupling constant). If the long-distance attraction due to pion exchange were sufficient to bind the 
pentaquark for any reasonable choice of short-distance dynamics (as happens in the combined limit) then one would 
have a robust prediction that heavy pentaquarks exist. Unfortunately, we find that this is not the case. 

Before proceeding it is useful to clarify a semantic point. Our discussion relies heavily on the large Nc limit of 
QCD; as Nc becomes large, the minimum number of quarks in a baryon containing a heavy antiquark is not 5, but 
rather Nc+2. Nonetheless, we still denote such states as "pentaquarks," to make the obvious connection to the Nc = 3 
world. 

This paper is organized as follows. In SeCn we provide a brief background on heavy pentaquarks. Section □ presents 
a rigorous argument for the existence of heavy pentaquarks in the combined large Nc and large niQ limits. In Sec. „ 
we discuss the symmetry structure of heavy pentaquarks in the combined limit, and in particular the fact that they 
fall into multiplets of SU(4) x 0(8) x SU(2). Then we explore in SeCnthe question of whether this qualitative result 
survives in the real world of Nc — 3 and finite mg by studying simple models based on a pion exchange between 
nucleons and heavy mesons. Finally, Sec. □ presents a brief discussion of the implications of this work and concludes. 

HEAVY PENTAQUARKS: BACKGROUND 

The experimental situation involving reports of heavy pentaquarks remains murky. The HI Collaboration at 
HERA has reported a narrow resonance Oc appearing in D*~p [{cd){uud)] and D*'^p [{cd){uud)] states produced 
in inelastic ep collisions, with a mass of 3099±3±5 MeV and a width of 12±3 MeV. We note that the Oc, even if 
it withstands further experimental scrutiny, is not the type of heavy pentaquark discussed in this paper, since it is 
a resonance unstable against strong decay. Moreover, subsequent evidence argues against its existence: The FOCUS 
Collaboration 21], using a method similar to that of HI but with greater statistics, finds no evidence for ©c. The 
experimental situation for heavy pentaquarks remains in a state as unsatisfactory as for their lighter cousins. 

On the theoretical side, much of the heavy pentaquark research to date has been performed in the context of 
different variants of the quark model ^13. 22, 23). Our purpose here is not to review this work in any detail, but to 
stress one of its key points: Heavy pentaquarks occur far more naturally than light pentaquarks in such models, simply 
because a heavy quark is drawn more closely than a lighter quark to the bottom of any potential well. At the time 
much of the theoretical analysis was performed, many researchers assumed that light pentaquarks were experimentally 
firmly established, and so such models seemed to make rather robust predictions of stable pentaquarks. Now that the 
existence of the light pentaquarks has become more questionable, the reliability of heavy pentaquark predictions can 
also be questioned. Nevertheless, the tendency of heavy pentaquarks to bind more tightly than light ones remains 
generically true, a simple fact that continues to play a crucial role in the analysis of this paper. 

Stewart, Wessling, and Wise jl^ also raise a critical issue in the context of a diquark type model, namely, whether 
heavy pentaquarks could prove stable against strong decays. They argue that negative-parity heavy pentaquarks 
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should have the lowest energy (in contrast to the positive-parity 0+ of the JafFe-Wilczek model '24|) since this involves 
s-wave interactions between the diquarks. They suggest that the additional attraction in such negative-parity states 
might be sufhcient to render the states stable against strong decays. In this paper we argue that pentaquarks do in 
fact exist, at least in the combined large Nc and large toq limits of QCD. 

Since the large Nc limit plays a critical role in our argument, it is useful to remark upon previous work on heavy 
pentaquarks as iVc^oo. References [i^H^li^ impose large Nc counting rules in the context of a quark picture as 
a way to implement large Nc QCD. Such a picture suggests a Hamiltonian and asymptotically stable eigenstates. 
However, generic excited baryons at large Nc are broad resonances with 0{Nc) widths and require an approach 
respecting their nature as poles occurring at complex values in scattering amplitudes. Two of this work's authors have 
developed just such a "scattering picture" jl^]. While obtainable through a generalization of the large Nc treatment 
for the stable ground-state band of baryons 0, the scattering approach naturally allows a proper treatment of 
resonant behavior such as large configuration mixing between resonances of identical quantum numbers . Even for 
pentaquarks of 0{Nc) widths, the scattering approach predicts multiplets degenerate in both mass and width [29l |. 
But this technology, while generally true, is not required in the current work; as we now show, the heavy pentaquarks 
discussed in this paper are stable against strong decay, at least in the combined formal limit Nc—>co, rriQ -^co. 



THE EXISTENCE OF HEAVY PENTAQUARKS 

We now show that heavy pentaquarks exist in the combined large Nc and large mg limits: They are stable against 
strong decay. We must first choose an appropriate parameter to describe the limiting procedure. Here, the natural 
choice is the A expansion, where 

A - 1/Nc , Aqcd/toq , (1) 

Aqcd is the hadronic scale, and mg is the mass of the heavy quark. We note that the natural expansion turns out 
to be in powers of A^/^ instead of A^ for a pure l/Nc expansion. 

Consider the states in the QCD Hilbert space that have energy less than + Mh + TOtt {Mh is the mass of the 
lightest hadron containing heavy antiquark Q), and have baryon number -1-1 and heavy quark Q number —1. These 
conditions exactly describe potentially narrow heavy pentaquarks Qq (assuming no symmetry forbids the one-pion 
decay). Now consider further states with energy less than Mn + Mh] any pentaquark state appearing here must be 
a bound state as no hadronic decay can occur. However, scattering states clearly occur between the nucleon and the 
heavy meson that have the appropriate quantum numbers and have low enough energies. Therefore states that can 
be labeled 9q exist. 

Yet is it possible to describe such a state as bound in some realistic potential? First note that momenta in the 
scattering states scale as A°. Therefore, since the N^ H reduced mass ^ scales as A~^, the kinetic energy scales as A^, 
which is much smaller than ttItt = 0(A*'). One may therefore construct an effective theory in which all scatterings with 
>2 final-state hadrons are integrated out. 

However, these states naively appear nonlocal, which would prevent the construction of a local potential. The range 
of the nonlocality scales as the inverse of momenta p associated with the smallest kinetic energy T one integrates 
out. In this case, T^m.,^. Therefore, the range scales like 1/p ~ (2/^TO7r)~^^^ ~ A^/^ as A^O: The nonlocality 
disappears. 

Next, one must ensure that the potential that binds the pentaquark does not vanish in the combined limit. From 
Witten's original Nc counting js^l, one finds that indeed V{r) ^ A", preventing its disappearance relative to the kinetic 
energy. Noting that the heavy quark coupling scales as Qs ^ Nc'^^'^ , the nucleon couphng is of order Qa/ f-n ^ Nc^^, 
and the pion propagator is of order ~ Nc, one combines these ingredients to find the desired A° scaling for the 
potential. 

We can now easily prove the existence of stable heavy pentaquarks. Having established the locality and scaling 
of the potential between heavy hadrons, we have successfully reduced a quantum field theory problem to one of 
nonrelativistic quantum mechanics. It is well understood in this context that a potential with an attractive region 
has an infinite number of bound states as /i — s- oo (see Appendix g for details). In the present case, ^ ^ A~^ ^ cxd, 
while y(f)~ A'^. Thus, proving the existence of heavy pentaquarks in the combined limit requires only that V{'r) is 
attractive in at least some region. Fortunately, we know the form of V{r) at large distances: It is given by a one-pion 
exchange potential (OPEP), because tt is the lightest hadron that can be exchanged between H and N. It is moreover 
known that, regardless of the relative signs of the coupling constants, attractive channels appear in the OPEP. Thus, 
V{r) necessarily has attractive regions, serving to bind the heavy pentaquark. 
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SYMMETRIES OF HEAVY PENTAQUARKS 

We now show that, in the combined large Nc and large mg limit, the pentaquark states form a multiplet of the 
group SU(4) X 0(8) x SU(2), which is an emergent symmetry of QCD. The SU(4) group is a spin-flavor symmetry 
of the Hght quarks similar to that in Refs. 0,13, 0|- The argument that SU(2)spinXSU(iV/)flavor combine 

to form a contracted SU(2A^/) is completely applicable to the case of heavy pentaquarks, where here we restrict to 
Nf = 2. 

The 0(8) group is the symmetry associated with the configuration of the heavy quark relative to the light degrees 
of freedom. For nonexotic baryons, the origin of this symmetry is explained in Ref. 18]. Since the reason for such a 
symmetry may not be so familiar, we provide further details here. Consider an attractive potential V{f) of the sort 
described in Sec. b Such a V^f) has a minimum, near which it can be approximated as harmonic. In the large Nc 
and large mg limits, the wave function is localized near this minimum, creating an emergent U(3) simple harmonic 
oscillator symmetry. This U(3) symmetry is generated by Tij = a|aj = 1,2,3), where aj is the annihilation 
operator in the jth coordinate direction. The generators satisfy U(3) commutation relations: 

[Tij,Tki] = SkjTu — SuTkj . (2) 

Additionally, as Nc oo the creation and annihilation operators also become generators of the emergent symmetry 
with the commutation relations 

[aj,Tki] = Skjtti , [a|, Tfei] = -Sual , [a^, a]] = % 1 , (3) 

where 1 is the identity operator. The sixteen generators {Tij, ai, a|,, 1} form the minimal spectrum- generating algebra 
for the U(3) harmonic oscillator. It is related to the U(4) algebra generated by Tij {i,j = 1,2,3,4) and satisfying 
commutation relations Eq. |(2J) by the limiting procedure 

aj ^ lim T4,j/R, a| = lim T4/R, 1= lim Ta/R^. (4) 

R — 'oc R — >oc R. — >oc 

Such a procedure is called a group contraction. Hence the group generated by {Tij, ai, aj., 1} is called a contracted 
U(4) group. 

The generating algebra of the contracted U(4) group can be expanded by including the operators Sij = OiOj and 
^ij ~ ^i*^] ~ ^' ^ith the following commutation relations: 

[Sij, Ski] = [Sik,ai] — 0, [Sij,al.] = SjkOi + SikOj, [Sij ,Tki] — SjkSu + SikSji, 

[Sij,Sli] = {SkiSij + SuSkj)'^ + SkiTj + SijTki + SuTkj + SkjTu , (5) 

while the commutation relations for Sjj can be obtained through Hermitian conjugation. This set of 28 generators 

{Sij,Sl^,Tij,au a^, 1} forms a closed operator algebra, which is a contracted 0(8). 

Reference [ig continues by showing that this emergent 0(8) is also an emergent symmetry of QCD. Extension to 
the present case is straightforward. The argument for the presence of the contracted 0(8) emergent symmetry relies 
on one's ability to approximate the bottom of V{r) as a harmonic oscillator potential. As we have seen, the large 
Nc and large mg limits ensure this feature by leading to fi—^oo. These conditions remain just as true for a heavy 
antiquark; thus the argument from |18j applies to heavy pentaquarks. 

The SU(2) is simply the symmetry of invariancc under spin rotations of the heavy quark: In the heavy quark limit, 
states with any alignment of the heavy quark spin are degenerate. 

BOUND STATES AND THE ONE-PION EXCHANGE POTENTIAL 

Now that we have shown stable heavy pentaquarks exist in the combined large Nc large toq limit, the critical 
question becomes whether they also occur in our Nc — 3 finite mg world. To our knowledge, this question cannot be 
answered in a model-independent way without solving QCD, and so we resort to models for enlightenment. 

We focus here on effective potential models based upon one-pion exchange at long distance. As discussed in Sec. „ 
such models are clearly useful not only because they represent physically correct phenomenology, but also guarantee 
stable pentaquarks in the combined limit. But we also note that the argument does not depend upon the particular 
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short-distance behavior of the effective potential. If the real world is sufficiently close to the combined-limit world for 
the argument to remain valid, all models of this sort must yield (multiple) stable pentaquarks. Note that the masses 
of the various pentaquark states can depend sensitively upon the details of the short-distance interaction, but their 
existence cannot. The question then becomes whether models of this type predict bound pentaquarks in a robust 
way, independent of the details of the short-distance physics. If so, one has a strong reason to believe that the states 
are, in fact, bound in nature. 

We construct a "realistic" potential that has the correct long-distance behavior (OPEP) and an ad hoc short- 
distance part constrained only by the natural scales of strong interaction physics. Our potential acts between a 
nucleon and a heavy meson (D or B). The nucleon-pion analogue is well understood; its interaction Lagrangian reads 

£Arjv. - -^^NT-j.^sNd^TT- , (6) 
/ttV 2 

where the axial coupling constant 5^ ~ 1.27, and the pion decay constant /jr^lSl MeV. 

However, the heavy meson-pion interaction is not as straightforward. We use a formalism similar to that outlined 
by Manohar and Wise 31] to encode the heavy quark symmetry. In the limit of Nq heavy quark flavors, QCD 
develops an emergent SU(2A^q) symmetry |32j . As a consequence of this symmetry, physical states do not depend on 
the spin of the heavy quark; thus the D{B) and D*{B*) mesons form a degenerate multiplet in the mc{mb) ^ 00 limit. 
The heavy meson-pion interactions can involve transformations between pseudoscalar and vector states {B <-> B* or 
D ^ D*). Using heavy quark symmetry, one combines them into a single field: 

il + ii) 

H^^—^[p;r~Pi,] , (7) 

where is the four- velocity, and the pseudoscalar and vector heavy meson fields are P and P*, respectively. This 
combination allows the interaction Lagrangian to be written in a manner similar to that of the nucleon interaction, 

£int = --^Tri7r"7^75i/a^^" . (8) 

/7rV2 

Of course, the pseudoscalar and vector mesons are not degenerate in the real world, due to l/mg corrections. The 
mass difference must be included in realistic models. 

Both the nucleon and heavy-meson interactions with the pion can expressed in terms of the spin and isospin of the 
particles: 

CNN.^'^^^{SN-d^'^)I% , (9) 



A„t = ^^^(s';.a^")/^ , (10) 

where Sn and In are the spin and isospin of the nucleon. Si is the spin of the light quark in H , and Ih is the isospin 
of the H field. Combining Eqs. © and Hl()|l . treating the nucleon and heavy meson in the static limit (i.e., ignoring 
recoil, which is suppressed in the combined limit) and Fourier transforming yields the OPEP in position space: 

14 (r) = lN-lH[2Si2VT{r)+ASN -SiV.ir)] 

= (/2 ll)[S,^VT{r) + {K^ -Si- S^)V.{r)] , (11) 



where the central part of the potential (r measured in units of l/m^r) is 

9A9H e" 
27r/2 3 



Vcir) = ^ V , (12) 



and the tensor part is 



T. / N 9AgH e 3 3 1 
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/ is the total isospin of the combined system, while K—Sn + Si, and 

Si2 = i[3{SN-f){Si-f)~SN-Si] . (14) 

It remains unknown whether and gn are of the same sign or of different signs, so the potential could have an 
additional overall negative sign. 

Clearly, the OPEP dominates the interaction at large r since tt is the lightest hadron. At shorter ranges the OPEP 
is no longer dominant and the effective potential is qualitatively different. The value of r at which the OPEP ceases to 
dominate the effective potential is presumably of order 1/Aqcd'^I fm, the characteristic range in strong interactions. 
Therefore, for distances less than some cutoff value tq ~ 1 fm, we use a purely phenomenological potential. Note that 
we do not simply add such a short-range potential to the OPEP at short distances, but entirely replace the OPEP 
by this new potential: The behavior of the tensor part of the OPEP at short ranges is unphysical and would 
completely dominate the potential if not removed. The short-distance potentials used are taken to be either (central) 
constants or quadratic functions, and their strengths are allowed to vary. If the logic of our underlying argument 
based upon the combined limit also holds for realistic toq values and Nc = 3, then the precise details of the potentials 
should be irrelevant to whether the pentaquark states bind. 

We use the OPEP of Eq. in a nonrelativistic Schrodinger equation and solve for bound states. Since the tensor 
term in the potential allows mixing between L states, L is not a good quantum number. However, S'12 commutes with 
the parity operator, making P a good quantum number. Therefore, states labeled by J, S (total spin S=Sq +K), 
and P are used as eigenstates. Treating states mixed under L requires a coupled-channel calculation; we obtain the 
coupled equations by including all possible states labeled by L and K that are consistent with a given set of J, S, 
and P. 

Lastly, since this potential is intended to be "realistic", in principle B-B* and D-D* mass differences can affect the 
results. Of course, these differences are 1/toq effects and vanish in the heavy quark limit. Since the principal reason 
for the model calculation is to test qualitatively whether we live in the regime of validity of the combined 1/Nc and 
l/niQ expansion, it makes sense to include this difference. However, in practice the effect of this mass difference is 
entirely repulsive, making the states are less likely to bind. Thus, if the states do not bind in the equal-mass case, 
they do not bind at all. Accordingly, we use equal masses and only investigate the effect of the mass splitting in cases 
where binding occurs. 

We attempt to make our model as realistic as possible, given the rather simple forms assumed for the short-distance 
potential. To this end, we choose for the heavy-meson coupling constant gn ~ ±0.59 (extracted from D* Dir 
decay [13; see below) and collect values for other observables in Appendix „ Table As an initial guess, we 
also constrain the parameters of the short-range potential such that this potential combined with a OPEP between 
nucleons gives the correct 2.2 MeV deuteron binding energy. This choice is not necessary, but it has the virtue of 
ensuring that the potential parameters are not completely unreasonable from the point of view of hadronic physics. 
We summarize the potentials in Tabled Ultimately, we vary many of the parameters in order to probe the robustness 
of the qualitative results. 

We then solve coupled differential equations using standard numerical methods. We seek bound-state solutions for 
all J = 5 and J = § states using both a constant and a quadratic form for the short-distance potential, for / = and 
1, and with either sign of gn relative to gA- Initially (as discussed above), we assume no mass splitting between 
the pseudoscalar and vector mesons. A complete set of tables of bound states thus obtained appears in Tables Hn and 
IIVI Here we focus on describing some key features of these results. 

For constant and quadratic potentials constrained by matching to the deuteron energy, bound states of the pen- 
taquark are quite sparse. No channel supports a bound state with a D meson. The B meson is able to bind weakly 
in the channels with negative parity, but only with / = 0. Binding in these states is relatively weak, around 1.3 MeV 
for the constant potential and around 3.9 MeV for the quadratic potential, and bindin g en ergies are consistently the 
same between these channels (Table iHll Cols. A and B). It should be noted that both [13 and our calculations have 
the negative parity states being more stable. The greater binding for the quadratic (versus the constant) potential is 
natural since it is significantly deeper. 

We also analyze the case in which the short-distance potential is simply set to zero. For this case, the OPEP does 
not bind a pentaquark for any channel. In order for this potential to bind without the aid of short-distance potential, 
gn would need to be raised to unreasonably high levels, near 1 (approximately double the extracted value), and in 
some cases larger than 2. When realistic mass differences between the vector and pseudoscalar mesons are introduced, 
binding becomes weaker. This mass splitting eliminates binding for all channels with either type of potential we 
consider. 

The heavy-meson coupling constant g^ used in our analysis is motivated by the results of a recent experiment by 
the CLEO Collaboration that measured |i33] the width of the D*^ D^n^ decay. The value of gn is extracted from 
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the width and found to be ±0.59 ± 0.07. The analogous decay process is energetically forbidden the in B sector, 
preventing a direct extraction; therefore, we exploited heavy quark symmetry and used the same value of gn for the 
B sector. Note, however, the uncertainty in the bottom sector due to possible l/mg corrections. Accordingly, we 
also investigated using a range of heavy-meson couplings and find the same qualitative results. 

These results depend upon the strength of the short-distance potential. Clearly, as these potentials become more 
strongly attractive, the states are more likely to bind. As the potential needed to bind deuterium may by anomalously 
small, a deeper constant potential was also considered. Table Hill Col. C and Table HVl Col. A show the results when 
the constant potential is decreased from the depth needed to bind deuterium, —62.79 MeV, to about 4 times as deep, 
—276 MeV. The deeper well both produces more bound states and causes previously unbound states to bind (In 
particular, the D meson can form a bound state in the deeper potential). 

The choice of OPEP cutoff at r = 1 fm is arbitrary. One does not expect the OPEP to be valid for r < 1 fm, but 
the effective cutoff might occur at somewhat larger r. Table HTll Col. D and Table Hvl Col. B present the binding of 
states with a cutoff of 1.5 fm (the potential depth is —62.79 MeV). The negative-parity states remain the only bound 
ones, but the binding is now stronger, and the D meson binds. These fluctuations in strength of binding indicate the 
importance of the short-distance physics to the heavy pentaquark formation. 

DISCUSSION 

Despite our general argument using the large Nc and large mg combined limit that the long-range OPEP is sufficient 
to bind pentaquarks, we find in our class of models that, if a heavy pentaquark binds at all due to one-pion exchange, 
it does so weakly in a few channels and depends in a nontrivial way upon the details of the short-range interaction. 
The main implication is obvious: In the real world, 1/A^c and l/mg corrections can be substantial. Indeed, they 
are large enough to render unreliable even qualitative predictions about heavy pentaquarks based upon the combined 
limit. 

Given this somewhat unhappy result, the most important question is whether or not heavy pentaquarks do in fact 
bind to form stable states under strong interactions, and if so, whether only very weakly-bound states occur, such as 
the ones seen here. Both of these questions remain open. We simply do not know enough about the short-distance 
part of the effective potential to provide a definitive answer. An optimistic view is that the short-distance interaction 
between the heavy meson and the nuclcon is likely to be more attractive than that between nuclcons, which has a 
strong repulsive core. This argument is particularly plausible if one views at least part of the repulsive core between 
nucleons to arise due to the Pauli principle between overlapping nucleon wave functions; this effect is greatly reduced 
in the interaction between a nucleon and a heavy meson. If it is true that the short-range effective potential between 
the heavy meson and the nucleon is significantly more attractive than the analogous nucleon-nucleon case, then it is 
quite likely that heavy pentaquarks form stable, tightly-bound states. 

Finally, we address the question of why the qualitative prediction of the combined large Nc and large mg limits 
is insufficient. At first sight this may seem surprising, since both the 1/Nc and l/mq expansions have proven to be 
predictive in many situations. One must remember, however, that the quality of a systematic expansion depends 
on coefficients as well as the expansion parameter, and the size of these coefficients depends on the observable 
being studied. If some observable has "unnaturally" large coefficients, then the expansion can easily fail unless the 
expansion parameter is extremely small. This view is echoed in The relevant question is whether one ought to 
expect "unnaturally" large corrections to the leading behavior. 

In retrospect, it is perhaps not so surprising that combined expansion is insufficient here. One can make an 
analogous argument, based entirely upon 1/Nc counting, that both the deuteron and the ^5*0 two-nucleon channel 
ought to be deeply bound and have a large number of bound states: Both the effective interaction between nucleons 
and the masses of the two nucleons grow as N^. However, as has been stressed elsewhere j^^, this argument fails 
for smaller values of Nc. Similarly, numerous doubly-heavy strongly-bound tetraquarks ought to exist in the heavy 
quark limit: The effective interaction between heavy mesons is independent of the heavy quark mass and scales as 
l/^Ncmq). However, as discussed in Ref. and based upon models similar to those studied here, it is questionable 
whether they are bound for finite mg . Evidently, the coefficients describing interactions between hadrons can in some 
qualitative way be sufficient to weaken significantly results one would naively expect directly from the 1 /Nc or 1 /mg 
expansions, yielding very large corrections to the leading-order results for real-world parameters. Why this is so is 
one of QCD's more intriguing mysteries. 

In conclusion, we showed that heavy pentaquarks must exist in combined large Nc and large mq limit, and that 
they form multiplets of SU(4) x 0(8) x SU(2). We constructed a one-pion exchange potential between a nucleon and a 
heavy meson, and solved coupled nonrelativistic Schrodinger equations, obtaining bound states. Some weakly-bound 



8 



states do exist in some models, but their existence depends upon unknown short-distance physics. The lack of binding 
emphasizes that the real world is too far from the idealized world of large Nc and large uiq to render the expansions 
robust for these observables. In order to deduce whether or not heavy pentaquarks exist requires a more complete 
understanding of the short-distance physics than is presently known. 
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BOUND STATES IN QUANTUM MECHANICS 

Consider a smoothly varying potential V{r) that vanishes as r^oo. If V{r) is nonsingular and has an attractive 
region, it must possess a minimum at some tq. In the neighborhood of ro the potential is approximately harmonic 
[i.e., V{r) ~ |(r-ro)2]. Therefore, if the wave function is for some reason localized near the minimum, then the 
system can be approximated as a harmonic oscillator. For large reduced mass the kinetic energy operator is small, 
and minimizing the wave function's curvature forces its localization near r = ro, as desired. The harmonic oscillator 
potential has an infinite number of bound states, separated by multiples of uj — ^kj [i. Thus we see that multiple 
bound states must exist for sufficiently large [i. If the potential is singular (but not more singular than l/r^, so that 
a ground state exists), the large size of [i localizes the wave function deep in the potential near the singularity, again 
allowing plenty of room for bound states. 

TABLES OF RESULTS 

This appendix focuses on our numerical results. Table lists the parameters used in the calculation. Table UTI 
summarizes the potentials that were used. Table Hill presents the energies of bound states for a B meson binding with 
a nucleon, while Table HVl presents the same for a D meson. 



Quantity Name 


Quantity Value 




1.27 




131 MeV 


9H 


± 0.59 




138 MeV 


ruN 


938.92 MeV 


niB 


5279 MeV 


mo 


1867 MeV 


As 


46 MeV 


Ab 


141 MeV 



TABLE I: Constants used in bound-state calculations for heavy pentaquarks. 
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K(^) 

Voir) 

Vrir) 
Vi(r) 
V2{r) 



-Il- lij)[S,2VT{r) + {K^ -Si- Sl)V.{r)] 
Vi{r) or V2{r) 

9A9H e""'-'' 2 
27r/2 3r " 



r > ro 
r < ro 



3 3 N 2 
2 2 "I h 1 



27r/| 6t 

VI) (Vo = -62.79 MeV or - 276 MeV) 
-252.659 + 541.321 r - 309.822 MeV 



fm 



(15) 



TABLE II: Potentials used in heavy pentaquark calculations. The labels are: total isospin 7, nucleon isospin 7jv, heavy meson 
isospin Ih, tensor force S12, tensor potential VrCr), nucleon spin Sjv, light quaxk in heavy meson spin Si, sum of nucleon spin 
and light quark spin K, central potential Vc(r). Numerical values are such that potentials are measured in MeV, distances in 
MeV~^, unless noted otherwise. Both Vi(r) and V-zi'r) arc central potentials. The parameters in V2{r) were fixed by making 
the potential differentiable at ro and bind deuterium with the appropriate energy. 



Channel 


I 


A 




B 


C 




D 




J 


S P 




+ 




+ 




+ 




+ 




1 
2 


1 
2 





1.30 


1.35 


3.89 1.92, 3.62 


139.38, 142.14 




14.49, 16.01 15.46, 16.15 






1 






0.35 


0.27 




139.38, 140.76 


15.32, 15.60 15.04, 15.46 


1 
2 


1 + 













14.9, 32.39 


4, 19.32, 46.5 










1 










12.72, 18.22, 26.91 


9.45 






1 
2 


3 
2 





1.30 


1.31 


3.89 


3.67 


140.76 


140.76 


15.87 


15.32 






1 








0.26 


140.76 


140.76 


15.04 


15.32 


1 
2 


2 ^ 













32.15 


3.35, 45.95 










1 










12.12, 27.19 


8.36, 22.08 






3 
2 


1 
2 





1.42 


1.31 


3.89 


3.67 


140.76 


140.76 


15.87 


15.32 






1 








0.26 


140.76 


140.76 


15.04 


15.32 


3 
2 


\ + 













15.32, 18.49, 32.43 


4.65 










1 










12.80 17.25, 17.66, 22.91 






3 
2 


3 
2 





1.42 


1.25 


3.89 


3.67 


140.76 


140.76 


15.87 


15.32 






1 








0.20 


140.76 


140.76 


15.04 


15.32 


3 
2 


2 ^ 













18.22, 32.29 












1 










4.18, 23.18 









TABLE III: B meson bound-state energies for each channel, where -I- and — refer to relative sign of qa and qh- All energies 
in McV. Column A: constant potential, Vq — —62.79 McV and ro = 1 fm; B: quadratic potential; C: constant potential, 
Vo = -276 MeV and ro = 1 fm; D: constant potential, Vo = -62.79 MeV and ro = 1.5 fm. 
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Channel 
JSP 


I 


+ 


A 


B 

+ 




1 1 

2 2 


— 




1 


113.99, 110.4 


114.82, 115.78 


7.36, 9.00 8 
8.40, 8.79 i 


i.45, 9.27 
i.l6, 8.63 


1 1 

2 2 


+ 




1 


2.91 


16 




1 3 

2 2 


— 




1 


117.3 

115.23 


116.2 
115.23 


9.00 
8.45 


8.45 
8.45 


1 3 

2 2 


+ 




1 


2.10 


15.87 






3 1 
2 2 


- 




1 


117.3 
115.37 


116.20 
115.78 


9.00 
8.45 


8.45 
8.45 


3 1 

2 2 


+ 




i 


2.91 




3 3 
2 2 






1 


117.3 
115.09 


116.20 
115.09 


9.00 
8.45 


8.45 
8.45 


3 3 
2 2 


+ 





2.53 












1 











TABLE IV: D meson bound-state energies for each channel, where + and — refer to the relative sign of qa and gn- AH 
energies in MeV. Column A: constant potential, Vo = —276 MeV and ro = 1 fm; B: constant potential, Vb = —62.79 MeV and 
ro = 1.5 fm. 



